Multi Shell Model for Majumdar-Papapetrou Spacetimes 
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Exact solutions to static and non-static Einstein-Maxwell equations in the presence of extremely 
charged dust embedded on thin shells are constructed. Singularities of multi-black hole Majumdar- 
Papapetrou and Kastor-Traschen solutions are removed by placing the matter on thin shells. Double 
spherical thin shell solution is given as an illustration and the matter densities on the shells are 
derived. 
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I. INTRODUCTION 

It was shown that the Majumdar-Papapetrou (MP) 
solution 0]- describes a spacetime possessing N ex- 
tremely charged black holes • Later, extremely charged 
static dust sources for MP spacetimes were considered 
and it was shown that solutions to Einstein-Maxwell 
equations with extremely charged dust restricted to thin 
shells exist. Moreover, solutions for various geometries 
including planar, spherical and cylindrical shells were ob- 
tained 5]. It was also shown that spherical shells can 
be used as sources for the Extreme Reissner-Nordstrom 
(ERN) spacetimes jBJ. The cosmological black- hole so- 
lutions which are time dependent generalizations of the 
MP solution were also studied 0- Q and generaliza- 
tions with discussions on non-radiative character of these 
spacetimes were discussed in 8] . 

In linear theories mass density of point particle distri- 
bution can be expressed as p(r) = ^ rrii6(r — ai) where 
m^'s are the masses of each particle, <5(r) is the Dirac 
delta function and a^s are the locations of the point par- 
ticles. This expression is in consistence with the Poisson's 
equation, 



A' 



V 2 V = -Anp with V = Y^ 



(1) 



i=l 



Such a point particle representation of the mass density is 
not valid in nonlinear theories. For instance, in a nonlin- 
ear theory like V 2 V + 4irpV 3 — 0, there exists no solution 
V for p = J2i Trii8(T — ai). We overcome the above prob- 
lem by replacing massive particles by spherical thin shells 
with centers at r = ai and radii r,;. Then the mass den- 
sity p can now be represented as p(r) = 53j=i POiS(Fi) 
where Fi(r) = (i = 1,2, .., AT) represents the positions 
of the shells and poi are the mass density of each thin 
shell. The space is divided into N + 1 spaces. In our 
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proposed method we shall solve such potential problems 
where the solutions are continuous everywhere. 

In this work, we will first review the solutions to 
Einstein-Maxwell equations in the presence of extremely 
charged static dusts restricted to thin shells that are 
given in Then we will generalize the methods in [{| 
to find solutions possessing extremely charged dusts re- 
stricted to multiple shells and obtain singularity free ver- 
sions of MP-solutions given in ||. We will also work 
out the double shell problem explicitly and show that 
the matter distributions on the shells are not uniform 
and that interior of the shells are not flat unlike the sin- 
gle shells in Finally we will work on cosmological 
MP solutions discussed in [||- 0- Q and develop similar 
multi shell models. 



II. SHELL MODELS FOR THE 
MAJUMDAR-PAPAPETRAU SPACETIMES 



Let M. be a four dimensional spacetime with the met- 



ric, 



ds 2 = -A 2 dt 2 + X 2 h il dx i dx j , 



(2) 



where hij is an Euclidean 3-metric and A is a function of 
spatial coordinates x l only. With this metric assumption, 
we consider a charged static dust source with the matter 
part of the energy- momentum tensor given by, 



rriM 



(3) 



where u u = — =- since the dust is static. The four- 
potential due to the static charged dust is in the form 
= Ao(x l )8^ l , so the current four- vector is given by, 



J» = a(x l )u^ = crix^XS ^ 



(4) 



where <j(x 1 ) is the charge density of the dust. Therefore 
the electromagnetic and the Maxwell energy-momentum 
tensors are given by, 



M. 



(5) 

~ -F a pF^g^). (6) 



2 



Using the Einstein tensor calculated from the metric (J2J 
and the field equations G M „ = 8tt(T*5+ m a^) with ©"© 
we get, 

K 

A' 



A 



K = ±1. 



(7) 



Then plugging the above equation back into the Einstein 
equations give, 



V 2 A + 47rpA 3 



0. 



(8) 



Therefore the Einstein-Maxwell equations reduce to a 
nonlinear type of Poisson's equation. To find the equa- 
tion satisfied by the charge density u{x l ) we use the 
Maxwell equations W V F» V = AirJ^ with |gj and © to 
get, 



V 2 A + 4ttkctA 3 = 0. 



Comparing (JSJ and © gives p — na which shows that 
the dust source is extremely charged. 

One can solve the reduced Einstein-Maxwell equations 
©-© for the extremely charged dust restricted to a thin 
shell as given in . Letting S be a regular surface in M 3 
defined by S = [(x, y, z) S M 3 ; F(x, y, z) = 0], the matter 
density can be written as, 



p(v) = p (v)6(F(v)), 



(10) 



where po is the matter distribution on the shell S. Gen- 
eralizing the method given in |f| we can set, 

A(r) = A (r) - Xx(r)Q{F) - A 2 (r)(l - 6(F)), (11) 

where 9(F) is the Heaviside step function. The above 
choice can be made since the spacetime M is now divided 
into two disjoint spacetimes M + and M~ by S where the 
metric functions A are Ao — Ai and Ao — A2 respectively. 
We also require that the following equations are satisfied, 



V 2 A = V 2 Ai = V 2 A 2 = 0, 
(Ai -X 2 )\s = 0. 



(12) 
(13) 



These follow from the fact that in M + and Mr the 
Laplace's equation is satisfied as can be seen from 
JHJ(Since p = Q in M + and M~) and the metric must be 
continuous across the shell. Inserting (|11|) into JSJl and 
using the conditions H12fl - I|13l) we S e tj 



Po = 



(VAi - VA 2 )-VF 



47TA 3 



(14) 



Let S be the sphere F = r — a = 0. One can choose 
A 2 = 0, Ai = A 3 (6>, - i) which clearly satisfies (|L2) - 
JI3J. Using (I3J gives, 



where A3 
1 



MM) 

satisfies, 



A, 



47ra 2 (A |s) 3: 



(15) 



d_f. Q d_ 

r 2 sin 89 V 89 



1 



8 2 



r 2 sin 2 9 8<p 



= 0. 



At this point we would like to remark that the above 
equation (fT*5f> can also be obtained by using the Israel 
junction conditions 0- However, when Israel junction 
conditions are used the surface energy-momentum tensor 
has the components, 



Stt — 



4™ 2 (A |s) 4 ' 



Sffff — Srhrh — 0. 



(16) 



So the matter density pi on S is given by using the fact 
that Stt = piu t u t = fr|s, 



A, 



PA 



4™ 2 (A |s) 2 



(17) 



As can be seen from <|17H in Israel method we get A 2 , in 



(9) the denominator instead of A 3 , as in which is due 



to a dimensional scaling. The total mass on the shell 
can be written as an integral of pod(F) over the 3 spatial 
dimensions which is equal to the integral of pi over the 
2 dimensional surface, that is, 



PoS(F) A d x = / pi\ z d z x, 



v 3 



(18) 



So 



So we get the same total mass on the shell in either 
method which means that there is no ambiguity. 

When the interior of the sphere S is chosen as a flat 
spacetime, we can set Ao = 1 and A3 = mo as constant 
which turns out to be the mass of the shell. So l|15|) 
becomes, 



Po 



(19) 



47TO 2 ' 

So the interior and exterior metric functions are given by, 

mo , toq 



A OM t 

A 7 77, 



1 - 
1. 



(20) 
(21) 



where too is the mass of the shell as can clearly be seen 
from Q19[l. This solution represents the ERN solution 
exterior to a spherical shell with flat interior. Letting 
r = (R — too) I j& where (3=1 — too /a with TOo^a assumed, 
the exterior metric can be written as, 



d.s- 



-/? 2 (1 



too, 
R ' 



'■dt 1 



to 



(l-^W + iiW, (22) 

it 



which is the conventional form of ERN metric (after a 
scaling). As we noted earlier, we aim to remove the sin- 
gularity of an ERN spacetime. Choosing too < a we 
see from (|20|l that \ out > is always satisfied. Com- 
paring 1|22JI. the relation r = (R — too)/ '(3 and 1)20(1 one 
can realize that the singularity of ERN spacetime is at 
X ut = (which corresponds to R — 0) and the horizon 
is at r — (which corresponds to R = too). This means 
that by restricting the matter on the shell with too < a 
we remove the both the singularity and the horizon of 
the ERN spacetime. In the case toq > a Q2Up. \ out = 



3 



is not excluded. The case mo = a represents the Levi- 
Civita-Bertotti-Robinson (LCBR) spacetime outside the 
shell and flat spacetime inside. By letting r = m^jR the 
exterior metric is obtained as, 



dR 2 + R 2 dn 2 }, 



(23) 



as the usual conformally flat LCBR metric. At this point 
we would like to note that in our analysis the limiting case 
a— >0 does not exist since in that limit the matter density 
HI U| l would become p(r) — po(r)5(r) and this choice is in- 
consistent with the Einstein-Maxwell equations JSJl"©- 
If we had the Newtonian theory, such a limit mathemat- 
ically would be consistent with the Poisson's equation. 
The non-existence of the a—>0 limit is in agreement with 
the results of |l0| . 

To remove the singularities of a MP spacetime pos- 
sessing N ERN black holes discussed in [3j , we place the 
matter source on thin shells. A MP spacetime possessing 
N ERN black holes has the metric function A given by, 



A = 1 



TV 

E rrij 
|r - ail 

3=1 1 



(24) 



where ai is the position of the i th ERN black hole. Such 
a spacetime contains N number of singularities. We can 
place the extremely charged dust on N spatially sepa- 
rated shells that do not intersect, instead of considering 
point sources. Generalizing the choice we made in 11U|) 
we can write the matter density of the spacetime as, 



TV 



(25) 



where poj is the matter distribution on the j th shell de- 
fined by Sj : Fj(r) = 0. We can also generalize the choice 
for A as, 

TV TV 

A = A - A e J] Q(Ej) - - e(f»], (26) 

3=1 3=1 

so that the metric function inside the j th shell is given 
by, 



Af = Ao-A, Vj = l,..,N. 



(27) 



Moreover the metric function exterior to all of the N 
shells is given by, 



\ ext \ \ 
A — Ao — A e . 

The functions Ao, Xj and A e satisfy the following, 



(28) 



V 2 A = V 2 Aj = V 2 A e = 0, (29) 



(A e -A,)| S3 =0 Vj = 1,..,tV, 



(30) 



where we again use the fact that the Laplace's equation 
is satisfied in source free regions and the metric across 



the shells must be continuous. Using H26() - (|29() and i|30|) 

we calculate V 2 A as, 



TV 



V 2 A = £(VA,- - VA e ) • VFMFi). 



(31) 



Using the fact that S(Fj)Q(F k ) 
given by, 



S(F 3 ) for k^j, p\ 3 is 



TV 



pX 3 



^^(Ao-A,) 3 ^). 

3=1 



Then inserting l|3*T|) and into © we get, 
(VA e - VAj) • VFj 



POj 



47t(A - Xj) 3 



(32) 



(33) 



As an illustration we consider the case of two spherical 
shells with radii r% and r 2 and the centers located at ai 
and a2. Our aim is to obtain the exterior MP solution 
with the metric function given in 124(1 so that the sin- 
gularities can be removed. The equations defining the 
spherical shells are given by F\ = |r — ai| — r\ =0 and 
F 2 = |r — a 2 | — r 2 =0. Making analogy with l(24|) and 
the single shell solution ((20(1 we choose, 



A R 



mi 



m 2 

T2 



mi 



m 2 



ail 



a 2 



(34) 



Appropriate choices for Ai and A 2 can be made by the 
following forms using and 



, to 2 
Ai — 



m 2 



Xo 



mi 

Tl 



|r - a 2 | 

mi 

l r - a i 



(35) 
(36) 



which in turn gives the full metric functions with the 
choice Aq = 1 as, 



Af 

A.r 



mi 


m 2 ^ 


n 


ti |r 


m 2 


m 2 


ri 


+ 1 r 

|r-a 2 | 


mi 


mi 


Tl 


+ |r-a x r 



mi 



m 2 



ai 



a 2 



(37) 
(38) 
(39) 



The above choices describe an exterior MP spacetime and 
ERN spacetimes inside each shell. Let ai = — ae z and 
a 2 = ae z so that the centers of the shells are on the z- 
axis with r± + r 2 < 2a. The surface matter distributions 
are obtained from l(3*3*) l with tV = 2 as, 



Pox 



P02 



mi 
Airr 2 



m 2 
A^ri, 



m 2 

T2 yfr\ 



mi 

Tl 



m 2 



4ari cos 9i + 4a 2 
nil 



\Jr\ + 4ar 2 cos 9 2 + 4a 2 



n -3 

(4P) 

-3 

(41) 



4 



where 6\ and 62 are spherical coordinates on Si and S2 
with their centers are taken as if origin. It can clearly 
be seen from (14011 and 1)41(1 that the existence of a sec- 
ond shell disturbs the uniform matter distribution on the 
other shell so that the system stays at equilibrium. The 
above equations can be obtained by Israel method as we 
discussed before, where one obtains the second power of 
the term in the denominator instead of the third power 
as in (|40|) and (|41(l which is due to dimensional scaling. 
One can obtain J23J) by letting 8=1 - &± - in {57)1 
so X ext becomes, 



X ext = 8 1 



mi 



m 2 



ai 



a 2 



= 8K 



(42) 



where fh\ = mi/ 8 and 7772 = m^/d. The MP space- 
time possessing two ERN black holes has two horizons 
each described as a surface enclosing the points where 
the black holes are located 0. These singularities are 
defined by the vanishing of X ext . So choosing the masses 
and the radii of the shells such that 777,1 /n + 777,2 /?"2 < 1 , 
we guarantee that X ext > from (|37|l . which means that 
by this choice we remove the singularities. Note that X™ 
in (|38[) and A|™ in (|39[l can never vanish in their domains 
of definition. The case m\/r\ + 7772 /r 2 > 1 is irrele- 
vant since we can not remove the singularities of the MP 
spacetime. The case mi/ri + 7772 /t^ — 1 should corre- 
spond to two mass generalization of the LCBR metric. 
Thus we obtained the double ERN black hole solutions 
with singularities removed. (Note that we obtained the 
solution with a re-scaling factor 8). At this point we 
remark that MP-solutions with multiple spherical dust 
shells as sources, can not have fiat interiors as we showed 
in 1|38|) and (|39|l . Moreover, the shells disturb each other 
and cause angular dependence of matter density on the 
shells. One can see from (|4(J|) and (|41|) that the second 
shell (S2) has maximum matter density at its north pole 
and the first shell(Si) has its maximum matter density 
at its south pole. Such a configuration puts the system 
in equilibrium. One can also calculate the masses of the 
shells by integrating l|4T)|) and (|4T|l on Si and S2 respec- 
tively as in l|18f) which gives the mass of the first shell 
as 7711 and the second shell as 7772 which is expected. We 
again note that the limits ri^O and ?"2— >0 do not exist 
since the matter density (|25|) is not consistent with the 
Einstein-Maxwell equations JSJ-©. This means that we 
can not have a source density of the form, 

p(r) = 777i(5(r — ai) + m 2 5(r - a 2 ). 

This is again in total agreement with |10| . 

One can also consider other solutions to the Laplace's 
equation for the metric function A to obtain various 
spacetimes. It was shown that such spacetimes possess 
naked singularities 0. One can remove these naked sin- 
gularities by placing the extremely charged source on thin 
shells as we did in this work for multiple ERN black hole 
solutions. The same procedure can also be applied to 
the stationary generalization of MP-spacetimes given by 



Israel- Wilson [llj and Perjes [T^ 
singularities of these spacetimes. 



to remove the naked 



III. SHELL MODELS FOR THE 
KASTOR-TRASCHEN SPACETIMES 

It was shown that time-dependent generalizations of 
MP spacetimes exist and solutions corresponding to N 
extremely charged co-moving black holes in a de-Sitter 
background were considered in 0- 0- The idea of 
removing the singularities of N co-moving black hole so- 
lutions by thin shells were considered briefly in |7fl for 
testing the cosmic censorship conjecture. At this point 
we would like to extend the discussion given in by the 
methods we discussed above. 

The cosmological MP solution which is a time- 
dependent generalization of the MP solution with metric 
PI in cosmological coordinates is given by, 



(j- 2 dt 2 + R(t) 2 U : 



h/q^j di*Xj .. 



(43) 



where R(t) is the scale factor and hij is an Euclidean met- 
ric. For a single extreme Reissner-N6rdstrom-de-Sittcr 
(ERNdS) black hole with its charge equal to its mass 
(Q = M) the metric function U in the metric (|43|l is 
given by, 



U = 1 



-Ht 



M 



(44) 



We can obtain the ERNdS solution in static coordinates 
by following the transformation given in [6( Given the 
metric ()43(l with U being arbitrary we consider a co- 
moving charged dust as source. As before, we assume the 
four-potential of the form = AqS^ and — ~|j<5° 
since the dust is co-moving (i.e. static in cosmological 
coordinates). Then the Einstein equations with positive 
cosmological constant become, 



(45) 



where M M „ is the Maxwell tensor given as in ©. The 
the Einstein equations l|4"5|) give, 



(-) 2 
K R J 

V 2 U 



u = 1 



A 

3' 

4:npR(t) 2 U 3 

1 



0. 



R(t) 



A(x), 



(46) 
(47) 
(48) 



where dot over R represents derivative with respect 
to t and A(x) is independent of t. The above equa- 
tions are obtained by separating the cosmological and 
electromagnetic-matter parts of the Einstein equations. 
One can note that l|46|l is just the Friedmann equation for 
a fiat cosmology. Solving l|4()l) yields R(t) = e Ht where 
H = ±y/Z/h. As discussed in @- 0- B negative H cor- 
responds to black hole spacetimes while positive H cor- 
responds to white hole spacetimes. We consider negative 



5 



H for the rest of the paper (H 
equations V ' V F^ V = iirJ^ give, 



\H\). The Maxwell 



V *U+—R(t) 2 U 3 = Q, 

K 



(49) 



with k = ±1 and a being the charge density of the co- 
moving dust. Comparing l|T7|) and l|4l3|) we get p = kg as 
before, so we conclude that the dust is extremely charged. 
Letting dr = R(t)dt we can write the metric l|43[) as, 



ds 2 = -U~ 2 dT 2 + U 2 h in dx l dx 3 . 



where 



U 
V 2 C/ 



Ht + X(x), 
inpU 3 = 0. 



(50) 

(51) 
(52) 



One can see from the Einstein equations 152fl that the 
product pU 3 must be time independent. This is clear 
from the fact V 2 [/ is time independent as can be seen 
from H51|) and to have the equation (|52|l satisfied for all 
times we must have pU 3 to be time independent. 

The above derivation of Einstein-Maxwell equations 
for a charged dust was given in 0. Our treatment of 
the shell model will be parallel to which we considered in 
the previous part of this work, that is different from the 
shells considered in . We start with a dust model and 
then take its limit to thin shells with the use of Dirac 
delta functions for the matter density. Since (|47|l and 
(|49|l are totally same with equations JSJ) and © we can 
consider similar multi-shell models with the matter den- 
sity given in (|25|l and the metric function U as in i|26|) 
(Clearly we just replace A with U). One difference is 
that the metric function U is time dependent via Ht. 
However, as can be seen from (|51|1 and (I52f) this time de- 
pendence does not affect Einstein-Maxwell equations. In 
source free regions we again have the Laplace's equation 
satisfied. The equation of each surface defining the shells 
are again given by Fj(r) — 0. Thus we follow the lines 
(EU to O to get the same dependence for the surface 
matter density. 

As an example, we again consider two spherical shells 
with radii r*i, r 2 and with centers located at ai, a2.Thc 
equations describing the shells are given by, 



F 2 



|r- ai| 
r — a 2 



n =o, 

r 2 = 0. 



(53) 



Due to linearity of Laplace's equation, we extend the 
single ERNdS solution in cosmological coordinates given 
in (|50|) to the double ERNdS case for the exterior metric 
function as, 



U' 



Ht 



m 2 
r 2 



nil 



m 2 



r - ai 



r-a 2 



(54) 



where we have chosen Uq = Ht in 128fl . Then 138|) and 

suggests, 



jjvr, 



Ht 



m 2 m 2 
r 2 |r 



a 2 



mi 

Ht 



nil 



ai 



(55) 
(56) 



for the metrics inside the first and the second shells. 
Clearly the choices H54l) - H55|) - (|56f) satisfy the required 
conditions (|29|l and 130JI for N = 2. We have mentioned 
above that pU 3 must be time independent which is sat- 
isfied in our multi-shell model since this product is given 
with the use of ifSKjl. (J2SJ) , JEJ and iJEIJ) as, 



PU 3 



p i(Ur) 3 S(Fi)+p 02 (Ur) 3 S(F 2 ) 



Since Fi t2 are time independent and poi, 2 Ul n 2 will be 
shown to be time independent below. 

Again choosing ai = — ae z and a 2 = ae z such that n+ 
r 2 < 2a the surface matter distributions are calculated 
from (|3H|l as, 



Poi = 



P02 



-i -3 



mi 


Ht - 


m 2 






r 2 


m 2 


Ht - 


m i 






ri 



rn 2 



r 2 yjr\ - Aari cos 0i + 4a 2 



mi 

ri ' ^rj+ Aar 2 cos 9 2 + 4a 2 



where the definitions of 0i and 9 2 are given as before. 
As can be seen from above, poi and po 2 are time de- 
pendent but when they are integrated on the surface of 
each shell (the term poi, 2 U 1 m 2 will appear in the integral 
which is time independent), they give constant masses 
for the shells as mi and m 2 which is in consistence with 
the result of Bonnor |]| . If we had a single shell, the 
interior of the shell would be de-Sitter. But for the case 
of double shells (and for N > 2 also) as can be seen from 
and (jnnjl the interior of the shells are ERNdS. Un- 
less this choice is made, the continuity condition of the 
metric across the shells will be violated. This fact was 
not realized in where the interior of the shells were 
considered as de-Sitter. 

To remove the singularities of the cosmological multi 
black hole solutions we again restrict the matter on thin 
shells. For that purpose, we assume that the thin shells 
are formed at r = with H < so from l|54|) one can see 
that JJ ext never vanishes for later times r>0. This means 
that the exterior metric is singularity free for all times. 
Moreover the interior of the shells described by (|55|) and 
(|56|l can also never vanish for the choice we made. Thus 
the whole spacetime becomes regular. The spacetime 
distances H-PQHg between two neighboring (space) points 
PandQis given by ||PQ|| 9 = |Z7|||P-Q|| where ||P-Q|| 
is the distance in R 3 , and U is given in l|51|) . Hence 
in exterior and interior spacetimes the space distances 
increase with constant (Hubble) speed. For this reason 
the shells inflate without collision, because the distances 
among them also increase. 

There is also another interesting possibility of con- 
structing regular spacetimes by considering H > for 
t<0 and H < for r>0. With these choices one can see 
that JJ ext , U[ n and U™ never vanishes and the spacetime 
metric is continuous at t — 0. In this model, the Hubble 
constant H changes sign in each universe separated by 
t = which causes a delta-function type of singularity 



6 



at this hypersurface. Spacetimes (interior and exterior 
spacetime regions) are contracting for t < and expand- 
ing for r > 0. 

Our method can be applied also to higher dimensional 
generalizations of the Majumdar-Papapetrau-Kastor- 
Traschen solutions 0] to study brane world creations 
and collisions 

IV. CONCLUSION 

We considered thin mass shell models for the 
Majumdar-Papapetrou and Kastor-Traschen spacetimes. 
In both cases we found exact solutions of the field 
equations describing gravitational fields of the extremely 
charged iV-number mass shells. Shell interiors are also 
curved spacetimes matching smoothly to the exteriors 
through (three dimensional) infinitely thin shells. Wc 
found the mass ( also the charge) densities on each shell. 
The solutions obtained this way are free of singularities. 



For the case of Kastor-Traschen spacetimes choosing the 
sign of the Hubble constant H properly we showed that 
shells are moving away from each other with constant 
velocity. 

Shell models in higher dimensional theories seem to be 
also very interesting. In order that our method to be ap- 
plicable to such theories one has to modify Einstein field 
equations with the inclusion of a dust matter as we did 
in this work where the mass is distributed on thin shells. 
In this case the shells are also higher dimensional. In the 
case of five dimensions for instance the thin shells are 
four dimensional spacetimes. Our work on this matter 
will be communicated elsewhere. 
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